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Abstract 



We perform a detailed study of the Earth matter effects on supernova neutrinos. 
The dependences of these effects on the properties of the original neutrino fluxes, on 
the trajectory of the neutrinos inside the Earth and on the oscillation parameters are 
described. We show that, for a large fraction (~ 60%) of the possible arrival times 
of the signal, the neutrino flux crosses a substantial amount of the matter of the 
Earth at least for one of the existing detectors. For oscillation parameters from the 
LMA solution of the solar neutrino problem the Earth matter effect consists in an 
oscillatory modulation of the v e and/or v e energy spectra. The relative deviation 
with respect to the undistorted spectra can be as large as 20 — 30% for E <; 20 
MeV and 70 - 100% for E Z 40 MeV. For parameters from the SMA and LOW 
solutions the effect is localized at low energies (E ^ 10 MeV) and is not larger than 
~ 10%. The Earth matter effects can be revealed (i) by the observation of oscillatory 
distortions of the energy spectra in a single experiment and (ii) by the comparison 
between the spectra at different detectors. For a supernova at distance D = 10 
Kpc, comparing the results of SuperKamiokande (SK), SNO and LVD experiments 
one can establish the effect at (2 — 3) a level, whereas larger statistical significance 
((4 — 5) a) is obtained if two experiments of SK-size or larger are available. Studies 
of the Earth matter effect will select or confirm the solution of the solar neutrino 
problem, probe the mixing U e 2, and identify the hierarchy of the neutrino mass 
spectrum. 
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1 Introduction 



Neutrinos from gravitational collapses of stars stop to oscillate in vacuum long time be- 
fore they reach the surface of the Earth. For mass squared differences implied by the 
experimental results on solar and atmospheric neutrinos, the supernova neutrinos arrive 
at Earth in mass eigenstates. The reason is either the loss of coherence between mass 
eigenstates on the way from the star to the Earth or the adiabatic conversion from flavour 
to mass eigenstates inside the star [|TJ. 

After a travel of thousands of years, the neutrinos can still be forced to oscillate again: 
it is enough to put on the way of the neutrinos a filter consisting of flavor non-symmetric 
matter, like the Moon or the Earth. The neutrinos will oscillate both in the filter and 
in vacuum after it. Indeed, the eigenstates of the Hamiltonian in matter do not coincide 
with the mass eigenstates z/j, i = 1,2, 3. So the latter turn out to be mixed and therefore 
will oscillate: z/j <-> Vj. As a consequence, oscillations in the flavour content of the flux 
are realized. 

The possibility of oscillations of supernova neutrinos in the matter of the Earth has 
been discussed long time ago It was marked that the effect of oscillations can be 
significant for values of parameters: Am 2 ~ 10 -6 — 6 ■ 10 -5 eV 2 and sin 2 29 > 2 ■ 10~ 2 . 
The effect is different for detectors with different trajectories of the neutrinos inside the 
Earth, and studying the oscillation effects in these detectors one can restore the direction 
to the supernova 0. 

The detection of the neutrino burst from SN1987A triggered a number of new studies 
of oscillations of supernova neutrinos and, among them, of oscillations inside the Earth. 
It was suggested that the difference in the Kamiokande-2 (K2) and 1MB spectra of events 
could be related to oscillations of v e in the matter of the Earth and to the different 
positions of the detectors at the time of arrival of the burst |J. For this mechanism to 
work one needs Am 2 ~ 10~ 5 eV 2 and large mixing of the electron neutrinos. 

The first two K2 events showed some directionality with respect to the supernova. 
The interpretation of these events as due to the scattering of electron neutrinos from 
the neutronization peak on electrons put strong bounds on the oscillation parameters, 
excluding a large part of the region which could be relevant for the solar neutrinos |4j]. 
In this region one expects the disappearance of the v e flux due to conversion v e — > v^, v T 
inside the star. It was marked however that oscillations inside the Earth can regenerate the 
v e flux, so that the bound is absent in the region of the large mixing angle (LMA) solution 
of the solar neutrino problem || || (see also the discussions in [J7[ ||). In this context, 
it was shown that the existence of the Earth matter effect depends on the conversion in 
the high density resonance associated with the large Am 2 which is responsible for the 
oscillations of atmospheric neutrinos || |K], |ll| . In particular, the adiabatic conversion in 
the higher resonance suppresses the Earth matter effects associated with the lower Am 2 
which governs solar neutrinos. 

At the cooling stage, when the fluxes of all neutrino species are produced, the conver- 
sion inside the star leads to partial or complete permutation of the v e and z/ M , v T spectra 
T\. This causes the appearance of an high energy tail in the v e spectrum which 
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contradicts the SN1987A observations Jl3| and therefore implies the exclusion of some 
range of the large mixing angles. However in the range of Am 2 ~ 10~ 5 eV 2 the bound is 
absent due to the Earth matter effect. 

Further studies of the Earth matter effect on supernova neutrinos have been performed 
in ref. [14 in connection with the role the supernova neutrinos can play in the reconstruc- 



tion of the neutrino mass spectrum. The three neutrino schemes both with normal and 
inverted mass hierarchy which explain the atmospheric and the solar neutrino data were 
considered. General formulas have been derived for neutrino and antineutrino fluxes at 
different detectors in presence of the Earth matter effect. The main qualitative features of 
the effect have been discussed and some examples of modification of the spectrum due to 
the Earth matter effect were given for schemes with LMA and SMA solutions of the solar 
neutrino problem. The role of the detection of the Earth matter effect in the identification 
of the neutrino mass scheme was clarified. In particular, it was shown that the very fact 
of the detection of the Earth matter effect in the neutrino and/or antineutrino channels 
will allow to establish the type of the mass hierarchy and to restrict the element U e s of 
the mixing matrix. 

In connection with the fact that the LMA gives the best global fit of the solar neutrino 

The Earth 



data, the interpretation of the SN1987A data has been revisited |15], [U], [17], [H 
matter effects on the antineutrino fluxes in the LMA range of the oscillation parameters 
have been studied in details (the two layer approximation of the Earth density profile 
has been used). The regions of the oscillation parameters have been found |16|] in which 
the Earth matter effects can explain the difference of the K2 and 1MB spectra. Such an 
interpretation also favors the normal mass hierarchy case or very small values of U e ^ for 
the inverted mass hierarchy [JT^, [16|, [L7J . 

where the expected 



Recently, the Earth matter effects were considered also in ref. |2( 
spectra of events at SuperKamiokande (SK) and SNO have been calculated in three neu- 
trino context with a two-layers approximation for the Earth profile. 



In this paper we perform a detailed study of the Earth matter effect on the supernova 
neutrino fluxes in the three neutrino schemes which explain the solar and atmospheric 
neutrino data. Using the realistic Earth density profile, we study the dependence of the 
regeneration effect on the oscillation parameters, on the trajectory of the neutrinos inside 
the Earth and on the properties of the original neutrino fluxes. The effects are calculated 
for both neutrinos and antineutrinos. 

We consider the possible arrival directions of the neutrino burst at existing detectors 
and estimate the probability to observe the Earth matter effect at least in two detectors. 
The possibility to identify the regeneration effect by existing and future experiments is 
discussed. 

The paper is organized as follows. In section we discuss the features of the neutrino 
fluxes originally produced in the star and the possible trajectories of the neutrinos inside 
the Earth. The general properties of the Earth matter effects are summarized in sect. 

a more specific discussion for LMA oscillation parameters is given in sect. [|. In sect. 
|5| we briefly discuss the regeneration effects for neutrino parameters from the LOW and 
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SMA solutions of the solar neutrino problem. The possibilities of observation of the Earth 
matter effects are studied in sect. ||. Discussion and conclusions follow in sections |7| and 



2 From the star to the Earth: trajectories and spec- 
tra 

2.1 Inside the star 

In this section we summarize the features of the neutrino fluxes as they are produced 
inside the star and the properties of the supernova that will be used in our calculations. 

A supernova is source of fluxes of neutrinos and antineutrinos of all the three flavours, 
e, /i, r. These fluxes, F° and F? (a = e,fj,,r), are characterized by the hierarchy of their 
average energies, 

(E e ) < (E g ) < (E,) , (1) 
and by the equality of fluxes of the non-electron neutrinos (which will be denoted as v x \. 

F* = Fi = F* = F« = F!. (2) 

In absence of neutrino mixing, 6 = 0, the neutrino flux at Earth is determined by the 
original flux produced in the star. If the decrease of the average energy and of the neutrino 
luminosity with time occurs over time scales larger than the duration of the burst the flux 
of the neutrinos of a given flavour, v a , can be described by a Fermi-Dirac spectrum as: 

L E 2 

E°(E, T a , L ai D) = eE/Ta - i , (3) 

where E is the energy of the neutrinos, L a is the total energy released in v a and T a is the 
temperature of the v a gas in the neutrinosphere. Here D represents the distance of the 
supernova from the Earth; typically D ~ 10 Kpc for a galactic supernova. The quantity 
E3 is given by F3 = 77r 4 /120 ~ 5.68. According to the hierarchy ([]]) the indicative values 
T e = 3.5 MeV, T s = 5 MeV and T x = 8 MeV will be taken as reference in our calculations; 
results will be presented also for other choices of the temperatures and their dependence 
on the specific values of T e , T E and T x will be studied. We consider equipartition of the 
energy between the various flavours, so that L a ~ Eb/Q, with Eb the binding energy 
emitted in the core collapse of the star: Eb — 3 • 10 53 ergs. 

In presence of neutrino mixing and masses the neutrinos undergo flavour conversion 
on their way from the production point in the star to the detector at Earth. Matter effects 
dominate the conversion inside the star, where a wide range of matter densities is met. 
The conversion effects depend on the distribution of matter in the star; the radial profile 

-to /10 Km\ 3 , , . 

,.(r) = 10"C(— ) g-cm" 3 , (4) 
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with C ~ 1 — 15 provides a good description of the matter distribution for p s 1 g ■ cm 3 
21], 0, |10|, p2[ . For p s ^ 1 g ■ cm -3 the exact shape of the profile depends on the details 



of the composition of the star. For the neutrino parameters we will consider the resonant 
transitions in the star occur at densities larger than 1 g • cm -3 , where the profile (^) 
applies; the value C = 4 will be used unless differently stated. 



2.2 Crossing the Earth 



If the neutrinos cross the Earth before detection, regeneration effects can take place due 
to the interaction with the matter of the Earth, analogously to what is expected for solar 



neutrinos 23 



What is the chance that one, two or even three detectors situated in different places 
of the Earth will detect the Earth matter effect on supernova neutrinos? 

Due to the short duration of the burst, and the spherical symmetry of the Earth, 
for a given detector the trajectories of neutrinos (and therefore the regeneration effect) 
can be completely described by the nadir angle 9 n of the supernova with respect to the 
detector: if cos 6 n > the detector is shielded by the Earth. The angle 9 n depends (i) on 
the location of the supernova in Galaxy, (ii) on the time t of the day at which the burst 
arrives at Earth and (iii) on the position of the detector itself. 

We first consider a supernova located in the galactic center (declination^ 5 S = —28.9°) 
and three detectors [P3| : LVD p5] , SNO |26| and SK ^7|. The positions of these detectors 
on the Earth are given in the Appendix A. 

The fig. [I] a) shows the dependence of cos^ on the time t for the three detectors. 
We fixed t = to be the time at which the star lies on the Greenwich meridian. The 
horizontal line at cos# n = 0.83 corresponds to the trajectory tangential to the core of the 
Earth (8 n = 33.2°), so that trajectories with cos# n < 0.83 are in the mantle of the Earth. 
For cos# n > 0.83 the trajectories cross both the mantle and the core. 

From the figures it appears that: 



1. For most of the arrival times the supernova is seen with substantially different nadir 
angles at the different detectors, so that one expects different Earth matter effects 
observed. 

2. At any time t the neutrino signal arrives at Earth, at least one detector is shielded by 
the Earth (cos# n > 0) and therefore will see the regeneration effect. Earth shielding 
is verified even for two detectors simultaneously for a large fraction of the times. 

3. At any possible arrival time t one of the detectors is not shielded by the Earth. So 
that, once the direction to the supernova is known, one can identify such a detector 
and use its data to reconstruct the neutrino energy spectrum without regeneration 
effect. 

1 We define S s as the the angle of the star with respect to the equatorial plane of the Earth. 
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4. For a substantial fraction of the times for one of the detectors the trajectory crosses 
the core of the Earth. 

In fig. [I] b), c), we show similar dependences of cos# n on the time t for other locations 
of the star in the galactic plane. Notice that in the case b) two detectors are not shielded 
by the Earth for most of the times and no one trajectory crosses the core. In contrast, 
for the position c) all the three detectors are always shielded by the Earth. In general 
for the detectors we are considering, which are placed in the northern hemisphere, the 
Earth shielding is substantial for a supernova located in the southern hemisphere, as it is 
the case for stars in the region of the galactic center. If a supernova event occurs in the 
northern hemisphere (5 S > 0), corresponding to some peripherical regions of the galactic 
disk, the Earth coverage of northern detectors will be scarce, or even null in the limit 
5 S = 90°. In this case southern detectors would be more promising for the observation of 
Earth matter effects. 

Clearly the determination of the position of the supernova is important for predic- 
tions and the experimental identification of the regeneration effect. The localization of 
the star can be done either by direct optical observations or by the experimental study 
of the neutrino scattering on electrons [^8|, which has substantial directionality. Trian- 
gulation techniques and neutron recoil methods have also been discussed [[28|, [29|. As we 
have mentioned in the Introduction, the study of Earth matter effects by high-statistics 
experiments can allow to reconstruct the direction to the star. 

3 Neutrino conversion in the star and in the Earth 

In this section we summarize the general properties of the Earth matter effect on supernova 
neutrinos. We will focus on the 3z/-schemes which explain the atmospheric and the solar 
neutrino data. 

3.1 Neutrino mass and mixing schemes 

We assume that the atmospheric neutrinos have the dominant mode of oscillations <-> v T 
with parameters [pPf ]: 

\m\ -m\\ = Ami = (1-5 - 4 ) ■ 10~ 3 eV 2 , sin 2 20^ > 0.88. (5) 

The solar neutrino data are explained either by vacuum oscillations (VO solution) or by 
one of the MSW solutions (LMA, SMA or LOW). The latter are based on the resonant 
conversion driven by the oscillation parameters 

\m\ — m\\ = Amq, sin 2 26*0 . (6) 

Moreover, we consider Am 2 atm 3> AmjL 

The electron flavor is distributed in the mass eigenstates v\ and v% with admixtures 
U e i ~ cos6 l 0, U e 2 ~ smO e . We will call i>\ and i>i the solar pair of states. 
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Three features of the neutrino schemes, which are important for the supernova neutrino 
conversion, are still unknown: 



1. The admixture U e z of the v e in the third eigenstate. Only an upper bound on this 
parameter is known from the CHOOZ and Palo Verde experiments [[H], |32fl : 



\U e3 \ 2 <0.02. (7) 

2. The type of mass hierarchy. In the case of normal mass hierarchy the solar pair of 
states is lighter than > m2,mi. In the case of inverted mass hierarchy the 
states of the solar pair are heavier than z/ 3 : 777,3 < m 2 ~ m i- 

3. The values of the solar parameters Airig, sin 2 26 , . Different solutions correspond 
to substantially different values of the oscillation parameters (see e.g. [ft3fl). 



The masses and mixings determine the pattern of level crossings in the star [|TJ] . There 
are two resonances (level crossings) in the schemes under consideration: 

• The high density (H) resonance, determined by the parameters Am 2 atm and U e3 . The 
conversion in the region of this resonance is described by the Landau-Zener type 
probability, Pjj, of transition between the mass eigenstates v 2 and is 3 . 

• The low density (L) resonance with parameters of the solar pair: Am^, sin 2 20 . 
We denote as Pl the probability of v 2 — * v\ transition associated to this resonance. 

Depending on the type of mass hierarchy and on the value of 9@, the resonances appear 
in different channels. There are four possibilities 



1. normal mass hierarchy and #0 < 7r/4: both the resonances are in the neutrino 
channel. 

2. normal mass hierarchy and 8 & > tt/A: the H resonance is in the neutrino channel, 
whereas the L resonance is in the antineutrino channel. This possibility is disfavored 
by the present data on solar neutrinos. 

3. inverted mass hierarchy and 9 & < ir/A: the H resonance is in the antineutrino 
channel, the L resonance is in the neutrino channel. 

4. inverted mass hierarchy and 9 & > 7r/4: both the resonances are in the antineutrino 
channel. 

These different schemes correspond to different conversion effects both inside the star 
and in the matter of the Earth. As it was shown in fT4]| , the Earth effects in the 3z/ 



context depend on (i) the type of mass hierarchy, (ii) the adiabaticity in the high density 
resonance, which is determined by \U e s\ and by the density profile of the star, (iii) the 
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oscillation parameters in the low resonance which are determined by the solution of the 
solar neutrino problem. 

In what follows we will consider the various possibilities in order. We will take oscil- 
lation parameters from one of the regions of the solutions of the solar neutrino problem, 
and also assume that the mixing parameter \U e ^\ is small, so that oscillations inside the 
Earth are reduced to 2v problem. 



3.2 Antineutrino channels 

Let us first consider the scheme with normal mass hierarchy. 



As discussed in sect. |3.1| , in this case there is no level crossing in the high resonance 
region in the antineutrino channel, so that the antineutrino flux at the detector does not 
depend on the jump probability Ph- We get: 



F g +(F°-F°)(l-2P i )(P le -|C/ el | 2 ) 



where 



F-e « Fl - (F° - F°)[(l - P L ) - (1 - 2P L )\U el \ 



(9) 



is the v e flux arriving at the surface of the Earth (without Earth matter effect) and the 
fluxes F® are defined in eq. (|3|). Here P\ e denotes the probability of V\ — > v e conversion 
inside the Earth and Pl is the jump probability in the L resonance. The forms @-© 
are the consequence of the approximate factorization of the dynamics and reduction of 



the three neutrino problem to an effective two neutrino conversion (see |14| for details). 



Let us consider the relative Earth effect expressed by the ratio: 

R EE £ " 6 ■ 



(10) 



From eqs. 



we find 



R = f(l- 2P L )f reg , 
where f is the ( "reduced" ) flux factor: 



e x 



P L + (l-2P L )\U el \ 2 +F° (l- j p L )-(l-2P L )|^ el | 2 



and f reg the regeneration factor: 



freg = (Pie ~ l^elP) • 



(12) 



(13) 



The three factors present in eq. ( |TT] ) describe the initial conditions (initial fluxes) and 
different stages of the evolution of the antineutrino state. Let us consider the properties 
of these factors in order. 
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1. The L-resonance factor, (1 — 2Pl), is close to the adiabatic value 1 (i.e. Pl — 0) 



especially if the resonance is in the neutrino channel So, in this case from eqs. 

O)-(O) we get the simplified expressions: 



R = ff,, 



reg j 



F||C/ el |2 + F0(l-|C/ el | 



(14) 
(15) 



2. The flux factor, eq. (|lq) , determines the sign and the size of the effect. Due to the 
hierarchy of energies, eq. (ED), a critical energy E c exists at which f = 0. Furthermore 
we have f > below the critical energy, E < E c , and f < for E > E c . For realistic 
temperatures of the neutrino fluxes (see sect. |2"TT| ) one gets: 

E c = (25 - 28) MeV . (16) 

At E ^> E c the flux factor ( jTj ) is dominated by the harder flux _F°, so that one finds the 
asymptotic behavior: 



f(£>>ft)= -(THbrw' (17) 



Similarly, at E <C £? c the flux F° dominates, giving the limit: 

1 



f(E <C £ c 



From eqs. (|T^)-(|T^) it follows that at very high, as well as at very low energies, the relative 
regeneration effect (0) becomes independent of the original fluxes. 



3. The regeneration factor, eq. (13), describes the propagation effect inside the Earth 



and is analogous to the regeneration factor which appears for solar neutrinos. Notice that 
f reg corresponds to genuine matter effect: it is zero in vacuum. 

The dynamics of propagation and properties of the regeneration factor (O) are differ- 
ent for oscillation parameters from different solutions of the solar neutrino problem. 

In what follows we perform numerical calculations of the Earth regeneration factor 
using a realistic density profile of the Earth We compare these results with results 
of the two layers approximation in the Appendix B. 

3.3 Neutrino channels 

If the hierarchy of the neutrino mass spectrum is normal, the H resonance is in the neutrino 
channel and the v e flux at the detector depends on P H JTJ 



' D - F e + (F e ° - F°)P H (1 - 2P L )(P 2e - \U e2 \ 2 ) , (19) 
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where the v e flux arriving at the surface of the Earth equals: 



F e ~ F° e - (F e ° - F°)[l - P H P L - Ph(1 ~ 2P L )\U e2 \ 2 } . (20) 
Here P 2e is the probability of the transition v 2 — > v e inside the Earth. 



From eqs. (|19|)-(p0|) one finds the relative Earth matter effect, R = (F® — F e )/F e , and 
the flux factor, r: 

R = rP H {l - 2P L )f reg , (21) 

r = ZJ. 12 l<y)\ 

F?P H [P L + (1 - 2P L )\U e2 \i} + F° [1 - P H P L - P H (1 - 2P L )\U e2 \>] ■ { ] 

The regeneration factor, f reg , is given by: 

freg = (Pe2 ~ \U e2 \ 2 ) = ~(P le - \U el \ 2 ) . (23) 

Let us comment on the features of the ratio R: 



1. From eq. (|2ll) it follows that if the adiabaticity in the high density (H) resonance inside 
the star is fulfilled, Ph — > 0, the Earth matter effect disappears. The reason is that in 
the adiabatic case the original electron neutrinos convert almost completely into and 
u T fluxes in the H resonance. Then the electron neutrinos detected at Earth result from 
the conversion of the original and v T fluxes. Since these fluxes are equal, eq. (@), no 
oscillation effect will be observed due to conversion in the low density resonance. 

The Earth matter effect is maximal in the limit of strong violation of the adiabaticity 
in the H-resonance: Pjj — > 1, when the dynamics is reduced to a two neutrino problem 
with oscillation parameters of the L resonance. 

The jump probability P H is determined by the density profile of the star and the os- 
cillation parameters \U e %\ 2 ~ tan 2 #13 and /S.m 2 atm . In fig. ^|we show the lines of equal Ph 
in the (Am 2 tm — tan 2 #13) — plane, together with the exclusion region from the CHOOZ 
experiment. We use the density profile (^); the error in Ph due to the uncertainty in 
the density profile is estimated to be within a factor of 2 Hli |. The figure shows that as 
\U e 3\ 2 decreases in the range allowed by the bound ([?[) the transition in the H resonance 
varies from perfectly adiabatic (Ph — 0), for \U e 3\ 2 ^ 5 ■ 10~ 4 , to strongly non-adiabatic 
(Ph — 1), for |t/ e 3| 2 ^ 10 -6 . The intervals of adiabaticity and strong adiabaticity vio- 
lation change only mildly as Am 2 tm varies in the presently allowed range. Notice that 
future atmospheric neutrino studies and the long base-line experiments will sharpen the 



allowed region of Am 2 



atm ' 



2. The low density resonance factor, (1 — 2P£), is zero if Pl = 1/2, which corresponds 
to a situation when the neutrino beam arriving at Earth consists in incoherent and equal 
fluxes of V\ and v 2 . In this case the effects of V\ v e and v 2 u e oscillations cancel each 
other. 

In fig. |3| we show the lines of Pl — 1/2, calculated with the density profile (f|) in the 
(tan 2 8 & — Am 2 & ) plane for different values of the neutrino energy. The lines cross the 
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allowed region of the SMA solutions for the energy interval E = 5 — 15 MeV. In fig. ||] we 
show also the lines Pl = 0.05 for two different energies: E = 5 MeV and E = 50 MeV. 
These lines determine the lower edge of the adiabaticity region in the (tan 2 # — Am^) 
plane. Notice that the LMA region is in the adiabaticity domain for all the relevant 
energies, whereas the SMA solution region is in the domain of significant adiabaticity 
violation. Depending on the details of the density profile of the star the LOW solution 
lies either in the adiabaticity region or in the region of partial adiabaticity breaking. 

A qualitative treatment does not depend on whether the low density resonance is in 
the neutrino or antineutrino channel (dark side of the parameter space). Quantitatively 
the results are different. 



3. The flux factor, r, eq. (f22|), changes sign at lower critical energy with respect to the 
case of antineutrinos, since the original v e spectrum is softer than the v e spectrum. We 
get: 

E c = (16 - 24) MeV . (24) 

Similarly to what was discussed for v e , the flux factor becomes independent of the original 
fluxes in the low and high energy limits: 

» - - i-p gft -fi!(i-2iw • (25) 

r(£<<£ ^ p H [^ + (i-2^ir (26) 



4. The Earth regeneration factor, eq. (|23D , depends on the mixing and mass squared 
difference of the solar pair and on the nadir angle 9 n . It will be described in detail in the 
following sections. 



3.4 Schemes with inverted mass hierarchy 

If the hierarchy of the mass spectrum is inverted the high density resonance is in the 
antineutrino channel (see sect. |3.1|) and the Earth matter effect for v e depends on the jump 
probability Ph- The expressions ©-(0) for the v e fluxes are immediately generalized to: 

F° = F- e + (F° - F°)P H {1 - 2P L )(P le - \U el \ 2 ) , (27) 
F- e « Fl - {Fl - F»)[l - P H P L - P H (1 - 2P L )\U el \ 2 ] , (28) 

in analogy with eqs. (|i"9|)-(pO|). 

Taking P L = (see sect. |3.2p , we find the relative deviation, R, and the reduced flux 
factor: 

R = rPnfreg , (29) 
• - , (30) 



F%P H \U el \* + F*(\-P H \U A 
11 



with the asymptotic limits: 



f(E >y ^ = - rw ■ (31) 
HE <K E ' } = ' (32) 

Clearly, the conversion of z/ e is independent of Ph- The expressions of the neutrino 
fluxes and F e can be obtained from eqs. fll9D -fl20D by the replacement Ph — > 1; they 
become analogous to eqs. With the same prescription, from eqs. (f2~T|)-(f2"2]) one 

gets the expressions of the ratios R and r. 

Summarizing the results of sections |3.2| - |3.4| we can say that the mass hierarchy and the 
adiabaticity in the H density resonance (and thus U e ^) determine the channel (y e or u e ) 
in which the Earth matter effects appear, which is: 

• both the v e and v e channels if the H resonance is strongly non-adiabatic, Ph = 1, 
regardless to the hierarchy. 

• the v e channel for adiabatic H resonance, Ph = 0, and normal hierarchy. 

• the u e channel for adiabatic H resonance, Ph = 0, and inverted hierarchy. 

The possibility of probing U e z and the mass hierarchy by the study of Earth effects on 
supernova neutrinos will be discussed in sect. [7]. 



4 The Earth matter effects for the LMA parameters 
4.1 Antineutrino channels 

Let us now study the features of the relative Earth matter effect, eq. fllij) , for u e with 
mixing and mass squared difference in the LMA region. 



The dynamics of the conversion inside the Earth is described by the regeneration factor 
f reg , eq. (|13|). For LMA parameters the Earth matter effect consists in an oscillatory 
modulation of the neutrino energy spectrum. 

The fig. |] shows the ratio R as a function of the neutrino energy for various values of 
9 n . For mantle crossing trajectories, 8 n > 33.2°, the effect is mainly due to the interplay of 
oscillations and adiabatic evolution. That is, oscillations in medium with slowly varying 
density. Small density jumps produce only rather weak effects. As a result, the features 
of the regeneration factor, f reg , are similar to what is predicted in the case of propagation 
in medium with constant density (see Appendix B). The factor is positive in the whole 
energy spectrum, so that the sign of the matter effect is determined by the flux factor 
(|1|): we have R > for E < E c and R < for E > E c . 
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The energy spectrum shows regular oscillations in energy with period 

AE « , (33) 

i, E) 



where the oscillation phase </> is determined by the integral 

Mn) = 2vr / . . n (34) 
Jo l m {n(x),E) 

over the neutrino trajectory. Here Z m is the (instantaneous) oscillation length in matter, 
and n(x) is the electron density along the trajectory. As follows from eqs. (|33T)- (p4|) the 
period of oscillations decreases with the nadir angle and increases with the energy. The 
dependence on the energy appears in AE, (see eq. fl33D ) explicitly, and implicitly via the 
oscillation length. 

As a result of adiabatic evolution, the depth of oscillations of the regeneration factor 
is determined by the electron number density at the surface of the Earth (see Appendix 
B),n°: 

D f « 2V2G F n° e -^ T sin 2 29° m . (35) 

Here (P m is the mixing angle of the solar pair in matter at the surface. 

The depth Df has a resonant dependence on the quantity x = 2E\V\/ Am^, with V 
being the matter potential (see the Appendix B for details). Both Df and l m increase as 
the system approaches the resonance; correspondingly, the period AE/E increases. For 
neutrinos propagating in the mantle and Am 2 Q = 5 ■ 10~ 5 eV 2 (which is used in the figure 
^) the resonance is realized at E — Er ~ 150 MeV. Thus the Earth effect is larger in the 
highest energy part of the spectrum. 

For core crossing trajectories the behavior of the Earth effect becomes irregular due 
to the interference of the oscillations in the core and in the mantle. The modulations in 
the energy spectra have smaller period both due to presence of large densities and larger 
length of the trajectory. Moreover, now the effect can change the sign both below and 
above the critical energy. That is, for some energies the Earth effect is negative at low 
energies and positive at high energies. 

As Am 2 Q decreases, as shown in fig. [5| the regeneration factor increases, since the 
resonance of the system is realized at lower energies. In particular, for Am 2 & = 2 • 10~ 5 
eV 2 , the resonance energy equals Er ~ 60 MeV. The period (in the energy scale) of the 
oscillatory modulation of the spectrum increases with the decrease of Am 2 Q . 

The change of the mixing parameter \U e i\ 2 ~ cos 2 (within the LMA region) in- 
fluences the regeneration factor rather weakly. However variations of \U e \\ 2 change the 
spectrum arriving at the surface of the Earth. According to the eq. (|9]), with the increase 
of \U e i\ 2 (i.e. decrease of the mixing 8 & ), the contribution of the hard component in the 
spectrum decreases: the composite spectrum becomes softer. 

Given the mixing U e \ the factor f, eq. (JT3J) , depends only on the original fluxes, F® and 
F®. The fig. ^| shows the relative Earth effect, R, as a function of the energy for different 
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values of f determined by different temperatures of the original v e spectrum. As it ap- 
pears in the figure, the critical energy E c decreases with T s , so that the region in which 
the flux factor suppresses the regeneration effect shifts to lower energies. For E <C E c and 
E ^> E c the depth of oscillations depends only very weakly on variations of T g , according 
to the limits fll8|) and (|17|) , which do not depend on temperatures. The values of T g and 
T x only affect the rapidity of the convergence to these limits: the convergence is faster 
for the larger difference T x — T g . This appears in fig. b), in which the flux factor f 
is plotted as a function of the antineutrino energy with the same parameters as in the 
fig. H For sin 2 2^0 = 0.75, as used in the figs. |-(7|, one finds f(E <C E c ) = 1.33 and 
f(E ^> E c ) = —4. Thus the Earth effect has stronger enhancement at high energies. 

If the hierarchy is inverted the Earth matter effect on u e is affected by the adiabaticity 
in the high density resonance. Such dependence is illustrated in fig. |8|, which shows the 
the relative effect R as a function of the energy for various values of Pjj. According to 
eq. ( |2D| ) the effect is proportional to Pjj and is maximal at Pjj = 1. The dependence of 
R on Pjj is transparent in the limits of high and low energies. Combining eqs. (p9|) and 
([31]) we find that at high energies R depends on Pjj as: 

R(E » E c ) = — P^yjreg , (36) 

1 — PH\Uel\ 

that is, R (x Pjj for Pjj <C 1. The weak dependence of R on Pjj in the softer part of the 
spectrum in fig. |] is explained by the fact that the dependence on Pjj cancels in the low 
energy limit, E -C E c , as can be seen from eqs. (^) and (|32|) . 



4.2 Neutrino channels 

Let us discuss the properties of the Earth regeneration effect in the v e channel. 

As it was shown in fig. |||, for LMA oscillation parameters the adiabaticity in the 
L-resonance inside the star is satisfied, so that Pj, = and eqs. (|21|)-(]22|) reduce to: 



R = rPjjf reg , (37) 

(38) 



F°P H \U e2 \2 + F°(l-P H \U e2 \ 



The regeneration factor f reg , eq. (|23|) , and therefore R, have similar dependence on 9 n and 
Am as in the case of antineutrinos. These dependences are illustrated in the figs. and 



10, where Pjj = 1 was taken. The oscillation length and the period of the modulations 



in the energy spectrum increase with the increase of the energy and the decrease of Am 
(fig. [T0[) . The depth of the oscillations of the regeneration factor f reg is larger than for 
antineutrinos since (if the L resonance is in the neutrino channel) matter enhances the v e 
mixing and suppresses the mixing of v e : 

sin 29 m {v) < sin 2 26 Q < sin 29 m (v) . (39) 
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The depth of oscillations has a resonant character (see Appendix B), increasing as the 
resonance energy is approached. According to eq. fl35|) the depth gets larger for smaller 
Ami ( fi g- 0)- 



The dependence of the Earth matter effect on the flux factor, r, eq. (|38|) , is illustrated 



in fig. O, where the ratio R is plotted for different values of the temperature T e , T x 



MeV and Ph = 1. According to sects. |37J and £373] the flux factor suppresses the Earth 
matter effect at energies close to the critical energy, E c , which is slightly lower than for v e : 
E c ~ 22 MeV for T e = 3.5 MeV. At high and low energies the asymptotic limits (p5|)-(p6|) 
are realized. This is shown in fig. ^ a), with the same values of the temperatures as in fig. 
TT| and Ph = 1. From eqs. (|25|)-(p6|) (with Pl = 0) and from the fig. ^ it follows that, 



in contrast to the case of v e , the Earth effect has stronger enhancement at low energy: 
for sin 2 2^ = 0.75 one gets r(E < E c ) = 4 and r(E > E c ) = -1.33. Notice that the 
convergence to these limits is faster than for v e due to larger difference of the v e and u x 
temperatures. 



According to eq. (0), the Earth matter effect is larger for larger Ph, i.e. for maximal 



adiabaticity breaking in the high density resonance inside the star. From eqs. (|37j) and 
the following asymptotic limit follows: 



R(E > E c 



P 



H 



\-P H \UA 2jre9 ' 



fr 



(40) 



similarly to eq. (p6|). The combination of eqs. (|37D and (^) implies that R becomes 
independent of Ph at E -C E c . These features of the dependence of R on P H are shown 
in fig. Q. 



As discussed in sect. d~4, the results for inverted hierarchy of the spectrum are obtained 
from the description given for normal hierarchy by the replacement Ph —> 1- Therefore 
the results shown in the figs. PHTT|, in which P H = 1 was used, apply to the case of inverted 
hierarchy. 



5 The case of oscillation parameters in the LOW and 
SMA regions 

5.1 LOW parameters 

For oscillation parameters in the region of the LOW solution of the solar neutrino problem 
the mass squared difference is smaller by at least two orders of magnitude with respect to 
the LMA solution: Am^ ~ 10~ 8 — 10~ 7 eV 2 . Therefore the resonance energy in the Earth 
is very small, Er < 1 MeV, and in the whole energy spectrum of supernova neutrinos the 
neutrino system is far above the resonance. Both the v e and v e mixings are suppressed 
in the matter of the Earth and the oscillation lengths approach the refraction length, 
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l m « Iq ~ 8000 Km. As a consequence, the regeneration inside the Earth has only weak 
dependence on the neutrino energy and no oscillatory distortions appear in the energy 
spectra at the detectors. In contrast, the regeneration effect depends strongly on the 
nadir angle 9 n : for E = 10 MeV the regeneration factor has a maximum, f reg ~ 0.04, for 
9 n ~ 25° (see e.g. @). 

The fig. [13] shows the relative Earth effects for v e and v e as functions of the energy 
for Am| = 10" 7 eV 2 , sin 2 2# = 0.9, 9 n = 25° and P H = 1. We have also considered 
perfect adiabaticity in the L resonance inside the star, Pl = 0. The figure shows that for 
E > 5 MeV the effect is not larger than 20% for u e and 10% for v e ; it decreases with the 
increase of energy with a <~ 1/E behavior. In the neutrino channel the relative deviation 
R, eq. (pip , is larger than in the antineutrino channel, especially in the very soft part of 
the spectrum, E ^ 10 MeV. This is explained (i) by the fact that the neutrino system 
approaches the resonance at low energies, and therefore the v e mixing is enhanced, and 
(ii) by the larger flux factor in the low energy limit, E E c (see fig. [7]). We get ~ 50% 
effect at E = 2 MeV and ~ 100% at E = 1 MeV. 



5.2 SMA parameters 

For oscillation parameters in the SMA region the fluxes at the detector, and F^, eqs. 
(H) and djjD, are substantially different with respect to the case of LMA parameters, due 
to differences in the L resonance factor, (1 — 2Pl), and in the regeneration factor. 

In what follows we describe the Earth effect for u e ; effects on v e conversion are ex- 
tremely small due to the suppression of the mixing in matter and will not be considered. 



Let us first discuss the factor (1 — 1P{). For AtUq and sin 2 29 & in the SMA region the 
conversion in the L resonance inside the star occurs in the adiabaticity breaking regime 
(fig. H). The jump probability Pl differs from significantly and, in particular, as shown 
in fig. |, Pl is close to 1/2, thus suppressing the matter effect (eq. (^Vj)). 



In fig. [14] we show the factor (1 — 2Pl) as a function of the neutrino energy for 
AttIq = 6 • 10~ 6 eV 2 and two values of sin 2 29 & in the SMA region. We used the profile 
(HI) with C = 4. For the largest possible mixing, sin 2 29 & ~ 5 • 10~ 3 , the factor is negative 
above 10 MeV; for the best fit point, sin 2 29 Q ~ 2.4 • 10" 3 , it is negative in the whole 
detectable part of the spectrum (E > 5 MeV). These results, however, strongly depend 
on the model of the star. For massive stars (M > 30M Q ) the density profile may have 
smaller gradient, so that the adiabaticity breaking is weaker. In this case the energy at 
which Pl = 1/2 is larger, and in a significant part of the spectrum the low resonance 
factor can be positive. Notice that in the high energy part of the spectrum the factor can 
be as large as 0.7 - 0.8 in absolute value. 



The regeneration factor, f reg , eq. (^), has a peculiar resonant behaviour in dependence 
on the energy and on the nadir angle 9 n . For mantle crossing trajectories, 9 n > 33.2°, the 
matter effect consists in resonantly enhanced oscillations, with oscillation length compa- 
rable or larger than the radius of the Earth. A peak appears in the u e spectrum at E ~ 15 
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MeV. 

For core crossing trajectories, 9 n < 33.2°, the regeneration factor exhibits a narrow 
peak at E ~ 7 MeV due to parametric resonance, with two smaller peaks at higher and 
lower energy due to MSW resonances in the mantle and in the core respectively. These 



features are shown in the fig. 15, which represents the relative deviation R as a function 



of the energy for Ph = 1, n = 0° and various values of the factor C in the density profile 
(H), corresponding to different values of Pl- With increase of C the jump probability Pl 
decreases, so that the region where (1 — 2Pl) > expands to higher energies. With the 
decrease of Pl the Earth matter effect changes from negative to positive and the size of 
the effect increases; it can be as large as ~ 30%. 

As in the case of the LMA solution, the effect decreases with Pjj and disappears for 
perfectly adiabatic transition in the H resonance, Pjj = 0. If the mass hierarchy is inverted 
the matter effect exists, in the neutrino channel, independently of the character of the H 
resonance. 



6 Observation of Earth effects 

Let us first summarize the results of the previous sections. 

The regeneration effects in the matter of the Earth produce a distortion of the neutrino 
energy spectrum. The character of the distortion depends (i) on the properties of the 
original neutrino and antineutrino fluxes, (ii) on the nadir angle 8 n and (iii) on the features 
of the neutrino masses and mixings, in particular on the solution of the solar neutrino 
problem. 

For oscillation parameters from the LMA solution the distortion of the spectrum has 
oscillatory character with larger oscillation depth and period in the high energy part of 
the spectrum. The effect exists in the neutrino or in the antineutrino channel or in both 
depending on the adiabaticity in the H resonance and on the mass hierarchy. 

For the LOW solution there is no oscillatory behaviour in the spectra. The effect is 
very small in the whole observable part of the energy spectrum and decreases with the 
increase of the energy. 

For parameters in the region of the SMA solution the Earth effect exists for neutri- 
nos only and consists in the appearance of one or three (for core crossing trajectories) 
resonance peaks (or dips) in the low energy part of the spectrum. 

In this section we discuss the possibility of detecting the Earth matter effects at present 
and future detectors. 

6.1 Detection of supernova neutrinos. Numbers of events 

As follows from the results of the previous sections, the observation of the Earth mat- 
ter effect requires: (i) separate detection of neutrinos of different flavours, (ii) separate 
detection of neutrinos and antineutrinos, (iii) the reconstruction of the neutrino energy 
spectrum. 
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In what follows we concentrate on events from charged current (CC) scattering on 
nucleons or nuclei at real-time detectors. These events better satisfy the requirements 
(i)-(iii); in particular CC processes have high sensitivity to the neutrino energy spectrum. 



We consider: 

1. The detection of v e at water Cerenkov detectors (SuperKamiokande and the outer 
volume of SNO) via the reaction: 

i? e + p —> e + + n . (41) 

Other CC reactions (e.g. the scattering of u e and v e on oxygen nuclei) have sub- 
stantially smaller cross section so that they contribute to the total number of events 
at few per cent level; they will not be considered further. 

2. Heavy water detectors (the inner volume of SNO experiment) with the detection 
reactions: 

u e + d^e + p + p, (42) 
u e + d-^e + + n + n, (43) 

which represent the dominant channel of CC detection. Events from the process 
(|43"|) will be distinguished by those from fl42|) if neutrons are efficiently detected in 
correlation with the positron. 

3. Liquid scintillator detectors (LVD), which are mostly sensitive to v e via the reaction 
(41) with only little sensitivity to absorption processes on carbon nuclei. 

Besides pure CC processes, the reactions 

Vi + e — > Vi + e , % = e, /i, r, (44) 
Vi + e — > Vi + e , (45) 

and the NC breakup of deuterium: 

z/j + d — > Vi + n + p , i = e, fi,r, (46) 
Ui + d-^^i + n + p, (47) 

allow to reconstruct the total neutrino flux; moreover, due to its good directionality, the 
scattering of v e on electrons is relevant to the location of the supernova. 

Radiochemical experiments could provide information on the total v e flux above a 
certain threshold. 

The number of CC events with lepton having the observed kinetic energy E e is given 

by 



ik =Nr £T dE ' ME - K)£(K) I dEFa{E) ^%F ' (48) 
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where E' e is the true energy of the electron (or positron), Nt is the number of target par- 
ticles in the fiducial volume and S represents the detection efficiency. Here cr(E' e , E)/dE' e 
is the differential cross section and R(E e , E' e ) is the energy resolution function, which can 
be described by a gaussian form: 



R(E e , E' e ) 



exp 



(Ee ~ K 

2A 2 



(49) 



The energy resolution A and the other parameters of the detectors (volume, efficiency, 
etc.) are summarized in the Appendix A. 

The energy spectrum (fS|) of the charged leptons reflects the spectrum of the neutrinos, 
with the following differences: 

• the energy dependence of the cross section^], a oc E 2 , substantially enhances the 
high energy part of the spectrum. 

• the integration over the neutrino energy and the convolution with the energy reso- 
lution function, eq. (f48[), lead to averaging out the fast modulations in low energy 
part of the spectrum (appearing for LMA oscillation parameters). Conversely, the 
large-period oscillations at high energies will appear in the lepton spectrum 



6.2 Identification of the Earth matter effects 

The Earth matter effects can be identified: 



1. at a single detector, by the observation of deviations of the energy spectrum with 
respect to what expected from conversion in the star only. 

2. by the comparison of energy spectra from different detectors. 

In the figs. |Tt3|4T9"| we show examples of the spectra expected at SK, SNO and LVD for 
oscillation parameters from the LMA solution, P# = 1 and various arrival times of the 
neutrino burst. We considered a supernova located in the direction of the galactic center 
(fig. [I] a)) at a distance D = 10 Kpc and releasing a total energy Eb = 3 ■ 10 53 ergs. The 
histograms represent the numbers of events from the reaction ( f4lf ) for SK (panels a)) and 
LVD (panels b)); the panels c) show the sum of the numbers of events from the reactions 
(f43|) and (fiTf ) at SNO. In d) we plot the numbers of events in the inner volume of SNO 
from the scattering (|42j). 

As can be seen in fig. |I] a), for t = 1 hour the neutrinos arriving at SK have core 
crossing trajectory (cos# n = 0.93). For SNO the trajectory crosses the mantle only and 
is rather superficial (cos# n = 0.10); the LVD detector is not shielded by the Earth. The 

2 Thc a oc E 2 dependence constitutes a good approximation at low energies; in the highest energy part 
of the supernova neutrino spectrum deviations due to weak magnetism and recoil effects are relevant, see 
e.g. |57l [}^, [39|. In our calculations we used the cross sections in ref. |39| for the scattering ( pd| ) and in 
ref. [[l0[ for the reactions (|^|)-(f43|). 
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corresponding spectra are shown in fig. [16] for Am?, = 5-10 5 eV 2 and sin 2 29 Q = 0.75 



and the same temperatures as in the figs. [| and |9| (see section |2.1| ). The spectrum of 



SK events exhibits deviations from the undistorted spectrum in some isolated bins, which 
correspond to the minima in the antineutrino spectrum. At SNO the Earth effect produces 
a narrowing of the spectrum. The difference with respect to the SK spectrum is attributed 
to the smaller oscillation phase (shorter trajectory in the Earth) of the neutrinos arriving 
at SNO. 

The comparison of the results of the three experiments can be performed in various 
ways. If the direction to the supernova is determined, the LVD spectrum is known to be 
free from regeneration effects. Therefore it can be used to predict the energy distributions 
at SK and SNO without Earth matter effects. Such predictions can be compared to the 
observations of SNO and SK. Due to the relatively small statistics of the LVD events, 
however, the accuracy of the reconstruction will not be high and the deviation from the 
undistorted spectrum (e.g. in the range 40 - 65 MeV) will not be larger than 2a. 

Higher statistical significance is obtained if data from a second large volume detector 
are available. Kilometer-scale neutrino telescopes, though primarily devoted to the study 
of high-energy neutrinos, are expected to be sensitive to supernova neutrinos. In particu- 
lar, for a supernova at distance D = 10 Kpc the ice Cerenkov detector of AMANDA jil) 
would observe more than 2 • 10 4 events from v e scattering on protons, eq. ( fill) H42 , fj3fl . 



Unfortunately, the presence of a relatively large background and the absence of sensitivity 
to the neutrino energy spectrum j|2], |3|, |4]| strongly restrict the potential of the study of 
this signal. We mark, however, that substantial upgrades of the experimental apparatus 
are possible Ji4 and the optimization of the detector for supernova neutrino observation 



would be of great interest. Besides the present neutrino telescopes, the detection of su- 



pernova neutrinos is among the goals of future large volume detectors, like UNO [45] and 



NUSL [55]. We find that the comparison of the energy spectra observed by SK and by 
another detector with comparable or larger statistics could establish the Earth matter 
effects at more than ~ 5 o level. 

Even larger regeneration effect can be realized due to specific features of the original 



neutrino spectra and of the solution of the solar neutrino problem. The figure [17] shows 
the same spectra as fig. [II] for Am = 3 • 10~ 5 eV 2 , sin 2 2^© = 0.9, T e = 3 MeV and 



T g = 4 MeV. In this case the regeneration effect is substantially larger (see also figs. [5] 
and |6|) and can be established by SK-LVD comparison with high (^ 3 a) statistical signif- 
icance. Notice also that the statistics of each experiment, and therefore the power of the 
comparison between different detectors, is higher for smaller distance D to the supernova 
and/or larger binding energy Eb- For instance for D = 3 Kpc and Eb — 4.5 ■ 10 53 ergs 
the statistics is ~ 17 times higher and the differences between the spectra of SK and LVD 
(unshielded) can be as large as (6 — 10) a. 

Besides the comparison of the spectra, more specific criteria of identification of the Earth 
effect can be elaborated if the location of the supernova and the solar neutrino oscillations 
parameters are known. For instance, for LMA parameters and rather superficial trajectory 
in the mantle the effect consists in a narrowing of the spectrum (see e.g. figs. [TBT-p 
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panels c)-d)). Thus the comparison of the widths of the spectra at different detectors 
may establish the Earth effect. 



As a further illustration, in figs. 18 and 19 we show the expected spectra for the same 



parameters as in fig. [16] but different arrival times of the signal (see fig. [TJ a) ) . For t — 8 
hours SNO is unshielded while SK and SNO have deep trajectories in the mantle. For 
t = 17 hours SNO has core crossing trajectory (cos^n ~ 1) and LVD observes effects of 
regeneration in the mantle only; SK is unshielded. 

For neutrino parameters from the SMA solution the Earth effects on the spectra are 
smaller than ~ 5 — 10%. This is explained by the fact that the narrow resonance peaks 
appear in the low energy part of the spectrum where the smoothing effect of the inte- 
grations (^) is stronger, thus suppressing the effect. Moreover, in the low energy region 
the detection cross sections and efficiencies are small. Therefore, regeneration effects with 
SMA parameters appear difficult to be observed. Similar conclusions hold for the case of 
LOW parameters, for which the effect is small and localized at low energies. 



As we have mentioned in sect. |3.4| , already the very fact of establishing the Earth 
matter effect in the neutrino and/or in antineutrino channel will have important impli- 
cations for the neutrino mass and flavor spectrum. For this it will be enough to study 
some integral effect of regeneration. Let us consider the following possibility. As we have 
discussed in sect. [|, for LMA parameters the regeneration effect is negative above the 
critical energy and moreover the relative size of the effect increases with E. The Earth 
matter effect at low energy is small mainly due to the dependence of the regeneration fac- 
tor f reg oc E (eq. (|3^)). Therefore to identify the regeneration effect one can compare the 
signals from the process fl4~l| ) in the various detectors at low, E < E s , and high, E > E, 



energies [47], where E s is some separation energy. Let us introduce the numbers of events 



Nl s i e - dE /jm , NmS r dE /2fM ( 5o) 



E th (IE, ' Je s (i /•;, 



e 



and the ratio 

n ^ir H - <«> 

In absence of Earth effects 1Z has the same value for every experiment provided that 
the detection efficiencies are independent of energy at E > E th , or in the particular 
case of equal efficiencies. Therefore, differences in the quantity ([|l]) are entirely due to 
regeneration effects: stronger effect corresponds to larger 7Z. In fig. [2(] we show the 
dependence of 7Z on the separation energy E s for different detectors; we have taken 
E t h — 20 MeV. The panels a) and b) of the figures refer to the situations illustrated 
in figs. [16] and [17] respectively. In agreement with the analyses of the spectra, in the 
latter case (panel b)) the effect is stronger: the deviation of the ratio 1Z for SK from 
the value 7Z in absence of regeneration, given by LVD, can be as large as (2 — 3) a. 
Larger deviation ((4 — 5) a) is realized if TZq is provided by an experiment with volume 
comparable or larger than the one of SK. 
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7 Discussion 



Within 2-3 years the solution of the solar neutrino problem can be identified by the 
results of the SNO, KamLAND and BOREXINO experiments. In particular, KamLAND 
will be able to establish the LMA solution and to measure AuIq and sin 2 26 & with 10 - 



20% accuracy [48], |5(| [51]. This will enormously sharpen the predictions for the Earth 
matter effects. 



The possibility exists that f/ e3 will be determined by MINOS |52| provided that its 



value is not too far from the present upper bound, eq. (|7]). In this case, U e3 is certainly in 
the adiabatic range. In a long perspective, a neutrino factory |53| will be able to cover the 
whole the range of U e ^ relevant for supernova neutrinos. Thus, either U e 3 will be measured 
or the upper bound on U e ^ will be so strong that the transition in the supernova will be 
certainly non-adiabatic. 

Let us consider possible implications of the supernova neutrino results depending on 
the solution of the solar neutrino problem. 

1). Suppose that the LMA solution will be identified with parameters close to the present 
best fit point (AttIq = (3 — 6) ■ 10~ 5 eV 2 ). As we have seen, in this case the regeneration 
effect can be observed. 

The features of the Earth matter effects depend on the value of U e ^ and on the type of 
mass hierarchy. For normal mass hierarchy and U e s in the adiabatic range (which appears 
as the most plausible scenario) we expect regeneration effects in the antineutrino channel 



and no effect in the neutrino channel (see sect. |3.4|) . In the supernova data further 
confirmations of such a possibility are (i) the absence of the neutronization peak in v e and 
appearance of the v^j v T neutronization peak, (ii) hard spectrum of v e during the cooling 
stage: (E e ) > (E s ). 

The relative size of the effect in v e channel, R, is determined by the regeneration 
factor and the flux factor, according to eq. (0). At high energies, when the flux factor 
reaches the asymptotic value (|TTD, eq. ([14]) gives: 

R(E » E c ) ~ -JZL- . (52) 
sm t/0 

That is, the effect is completely predicted in terms of solar oscillation parameters. 

In the case of inverted mass hierarchy the Earth matter effect should be observed in the 
neutrino channel and no effect is expected in the antineutrino channel if \U e ^\ 2 > 10~ 5 
(see fig. 0). This possibility will be confirmed by the observation of the i/ e -neutronization 
peak and of an hard spectrum of the v e during the cooling stage. The relative size of the 
Earth matter effect, R, is given by eq. (l37f) (with the replacement Ph — > 1) with the high 
energy asymptotic (see eq. (fRj)) 
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In the limit | 1 2 <C 10~ 5 the high density resonance is inoperative, so that the result 
is insensitive to the mass hierarchy. Oscillations appear in both the neutrino and an- 
tineutrino channels, and, at high energies, they are determined by the solar oscillation 
parameters (eqs. (|52| ) and ([53])). The ratio of the relative effects at high energies equals: 

f^tan 2 £ ^. (54) 

Jreg 

So, possible checks of this equality would be the confirmation of the neutrino scheme with 
very small 

If U e z is in the intermediate region: \U e 3\ 2 ~ 10~ 6 — 10 — 5 the situation is more complicated. 
One expects to observe oscillations both in the neutrino and antineutrino channels; the 
regeneration effect depends on the mass hierarchy and on the specific value of U e %. In the 
case of normal mass hierarchy, the relative effect in the neutrino channel is proportional to 
Ph = PH(U e 3) (eq. (137|) ) and at high energies, when the flux factor reaches the asymptotic 
value we get, (eq. fl40D ): 

R{E>>E ^- i-p"^eJ *"- (55) 

Since f reg is determined by the solar parameters, by measuring the relative deviation R 
we can determine the value of Ph via the eq. (|55| ) and therefore get information about 
U e3 . 

For antineutrinos with inverted hierarchy we find the high energy asymptotic (see eq. 

R(E » E c ) ~ Ph f reg . (56) 

1 — Ph cos^ q 



In practice, the observation of the Earth matter effect in the u e channel and absence of 
the effect in v e channel will testify for normal mass hierarchy and \U e ^ > 10~ 5 . In the 
opposite situation, effect in the v e channel and absence of the effect in v e channel, the 
inverted hierarchy will be identified with |£7 e 3| 2 > 10~ 5 . However the present experiments 
have lower sensitivity to v e fluxes with respect to the fluxes of v e , so that it may be 
difficult to establish "zero" regeneration effect with high enough accuracy. 

If the Earth matter effect is observed in both channels, one should compare the size 
of the effect with that predicted in the absence of the high resonance in a given channel. 
Thus, if the observed signal in the neutrino channel is smaller than what is predicted in 
the assumption of Ph = 1, whereas in the antineutrino channel prediction and observation 
coincide, we will conclude that the hierarchy is normal and the ratio of the observed to 
predicted signals in the neutrino channel can give the value of Ph- The opposite case of 



23 



coincidence of the predicted and observed signals in the neutrino channel and suppressed 
observed signal in the antineutrino channel will testify for the inverted mass hierarchy. 

Besides the probing of the neutrino mass spectrum and mixing, a study of the prop- 
erties of the original neutrino fluxes can be done with Earth matter effects. In principle, 
a detailed study of the observed energy spectra will allow to reconstruct the flux factor 
as well as to determine the critical energy E c . 

2). Suppose that the future solar neutrino experiments will identify the SMA solution. 
In this case, the Earth matter effect is expected in the neutrino channel only, and only 
if the high density resonance is inoperative. This requires the inverted mass hierarchy or 
very small U e 3 in the case of normal mass hierarchy. As discussed in sect. |5.2| , the effect 
is small and difficult to be observed even in the most favorable situations. 

For the rather plausible case of the normal mass hierarchy and \U e s\ 2 > 10" 5 no Earth 
matter effect should be seen. 

The observation of the Earth matter effect in the SMA case will allow to conclude 
that the mass hierarchy is inverted or the hierarchy is normal but U e s is very small: 
| C/ e 3 1 2 <C 1CT 5 . In principle, the intermediate case \U e s\ 2 ~ 1CT 5 can be identified if the 
observed signal will be smaller that the expected one for Ph = 1. Notice that from eqs. 
and (|22|) one gets the low energy asymptotic for the flux factor: 

R(E « E c ) = P H (1 - 2P L )f reg -^— = f reg (l - 2P L )i- , (57) 

-TfT-TL U 

which does not depend on P H . So it will be difficult to disentangle the effect of Ph from 
uncertainties in the original neutrino fluxes. 

The situation can be much more complicated if a sterile neutrino exists. This can 
be clarified in 2 - 3 years: the MiniBooNE experiment and further searches for sterile 
neutrinos in the solar and atmospheric neutrino experiments will allow to establish the 
existence of sterile neutrinos. 

Negative results of the searches will strongly favor the 3z/ schemes discussed in this 
paper. Still some uncertainty will remain: sterile neutrinos, unrelated to the LSND result, 
may exist and weakly mix with active neutrinos. Even a very small mixing (unobservable 
by other means) of sterile states with masses in the wide range from sub eV up to 10 keV 
can strongly modify the properties of the neutrino burst. 

8 Conclusions 

1). There is a big chance that at least one of the existing detectors (SK, SNO, LVD) will 
be shielded by the Earth at the moment of arrival of a supernova neutrino burst, so that 
the Earth matter effect on the neutrino flux will be observed. For supernova in a region 
close to the galactic center the most plausible configuration is that for two detectors the 
trajectories of the neutrino burst cross the Earth, whereas the third detector is unshielded. 
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We found that the detectors considered can register Earth matter effects (cos# n > 0) 
for a significant fraction 60%) of the possible arrival times of the signal and core effect 
for ~ 20% of the times. These fractions may be even larger depending on the specific 
location of the star in the galaxy. 

The comparison of the signals from different detectors allows to identify the Earth 
matter effect and to get information on the neutrino mass spectrum, substantially reduc- 
ing the uncertainties related to the model of the star and to the original neutrino fluxes. 

2) . We studied the effects of the matter of the Earth on supernova neutrinos in the frame- 
work of three flavours with either normal or inverted hierarchy. 

3) . The strongest regeneration effect is expected for oscillation parameters in the region of 
the LMA solution of the solar neutrino problem, especially for the lowest values of Am| 
in this region: Am@ = (2 — 5) • 10~ 5 eV 2 . 

In the v e — channel the effect exists in the scheme with normal mass hierarchy (ordering 
of the states) or in the scheme with inverted mass hierarchy for \U e ^\ 2 < 10 -5 , when 
the conversion in the high density resonance is non-adiabatic. The effect consists in 
an oscillatory modulation of the energy spectra and is negative (except in small energy 
intervals for core crossing trajectories) above the critical energy E c ~ 25 MeV, thus 
suppressing the signal. The relative size of the effect increases with energy and at E ~ 
60 — 70 MeV it can reach 50 — 70%. The period of modulation increases with energy 
and above E ~ 40 MeV no averaging occurs in the energy spectrum of events. At low 
energies the effect is small and the modulations are averaged out. Thus, in the LMA case 
the most sensitive region to the Earth matter effect is above E ~ 40 MeV . 

The oscillatory picture (position of minima and maxima) is very sensitive to AtUq. 
For trajectories in the mantle only the modulation of the spectrum is rather regular. For 
the core-crossing trajectories the structures become narrower and irregular. 

The Earth matter effect decreases with the increase of Am 2 & and for Am| > 10~ 4 eV 2 
it will be difficult to be observed. 

4) . For the v e — channel substantial Earth matter effect is expected in the case of normal 
mass hierarchy, provided that |^7 e3 | 2 < 10~ 5 , or in the case of inverted mass hierarchy. 
The effect has oscillatory character, similarly to the v e — case, and can be as large as 100 
% at E > 50 - 60 MeV. 

Thus, as in the case of antineutrinos, one should search for an oscillatory modulation 
of the signal in the energy range E > 40 MeV. 

The size of the effect depends on the properties of the original neutrino fluxes via the 
flux factor. This dependence, however, disappears in the high-energy limit, E > 50 MeV. 

5) . For oscillation parameters from the SMA solution the effect appears in the v e — chan- 
nel only, and only if the H-resonance is inoperative (very small C/ e 3 or inverted mass 
hierarchy). The relative effect can be as large as 30% for core crossing trajectories due 
to the parametric enhancement of oscillations. The effect is localized in the low energy 
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part of the spectrum: E = 5 — 10 MeV. It can be further suppressed by the L-resonance 
factor (1 — 2Pl) which in turn strongly depends on the density profile of the star in the 
outer region. The effect will be strongly smoothed in the spectrum of observed events 
by integrations over the neutrino energy and the true energy of electron. Practically no 
distortion is seen and the effect is reduced to an increase of the number of events in the 
region of the peak. Effectively this will make the spectrum narrower. One can probably 
identify this effect by comparing the signals from two high statistics experiments. 

6) . In the case of LOW solution the Earth matter effect is significant at very low energies 
with a smooth 1/E behaviour. If the H-resonance is inoperative, in the v e channel the 
effect (excess of flux) can reach 100% at E — 1 MeV and 20% in the v e channel (the 
difference is mainly due to the difference in the flux factors). At E = 5 MeV, the effect 
in the flux is below (10 — 20)% . The nadir angle dependence is determined by the length 
of the trajectory: the oscillation length almost equals the refraction length and depends 
only weakly on the neutrino energy. For the core crossing trajectories the effects can be 
enhanced due to parametric effects. 

7) . The identification of the Earth matter effect is possible in a single detector (in the 
LMA case) by the observation of the oscillatory modulation of the energy spectrum in 
the high energy part: E > 40 MeV. Another method consists in the comparison of signals 
from two (or several) different detectors. If one of the detectors is unshielded by the Earth 
its result can be used to reconstruct the spectrum of the neutrinos arriving at Earth and 
make predictions of the signal expected in the other detectors in absence of matter effect. 
For a supernova at distance D = 10 Kpc and with energy release Eb = 3 • 10 53 ergs we 
estimated that the Earth matter effects can be established at (2 — 3)<r level by comparison 
between SK, SNO and LVD results, and at (4 — 5)<x by comparison between the spectra 
from two large volume detectors (of SK size or larger). 

Another method to identify the Earth matter effect is to study the ratio of numbers 
of events in the high and in the low parts of the spectrum in different detectors. 

8) . Studies of the Earth matter effect will allow to establish or most probably to confirm 
the solution of the solar neutrino problem, to get information about U e s and to identify 
the type of hierarchy of the neutrino mass spectrum. 
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Appendix A. Parameters of detectors 

We summarize here the characteristics of the SK, SNO and LVD detectors that have been 
used in the analysis of section |6.2| . 



For each experiment we consider: 



1. The position of the detector on the Earth, which is relevant for determining the trajec- 
tory of the observed neutrinos inside the Earth (see sect. |2.2| ). The locations of the three 
experiments are given in Table [I] in terms of northern latitude, 5, and eastern longitude a. 

2. The material which the detector is made of and its fiducial mass. These quantities are 
quoted in Table [l| 



3. The detection efficiency S(E' e ) (see eq. 

The SNO efficiency is high, so that the shape of the energy spectrum and the total 
number of events are determined by the detection cross section [|53 |. Therefore, we have 



taken 8 = 1 in eq. ( 48 ) 



The efficiency of the LVD detector has been provided by the dedicated collaboration 



| 55fl . For E t h < 10 MeV, it is given by the gaussian integral function: 



£(E m ,E th ) = — L= / e-^dy, (58) 



/2tt 

where x = (E m — E t h)/^th an d E m is the measured energy of the lepton. The values of 
the parameters E th and a t h are: E th = 4 MeV and a t h = 0.74 MeV for the core counters 
(mass M = 0.57 ktons); E t h = 7 MeV and ath = 1-1 MeV for the external counters (mass 
M = 0.43 ktons). 



For SK we adopted the efficiency published for the Kamiokande-2 experiment [56 



4. The energy resolution A, which appears in the resolution function (|4"9|). We followed 
ref. |57| for SK and SNO experiments, and refs. [58], |55| for LVD. The parameter A 



depends on energy as follows: 



MeV MeV V MeV 
The values of the coefficients A and B are given in Table 



A „ E e / E e 

A—\- + B\—\-. (59) 



Appendix B. The regeneration factor: step- like and 
realistic Earth profile 

The analytical expressions for the regeneration factors (p3|) and (|i~3"D can be obtained in 
the two layers approximation of the Earth density profile. In this approximation the man- 
tle and the core of the Earth are considered as layers with constant densities. Therefore 
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detector 


5 


a 


material 


mass (ktons) 


A 


B 


SK 


36°21' 


137°18' 


H 2 


32 





0.5 


SNO 


46°30' 


-sroi' 


H 2 


1.4 





0.35 








D 2 


1.0 






LVD 


42°25' 


13°31' 


(CH 2 )io 


1.0 


0.07 


0.23 



Table 1: Summary of the characteristics of the detectors we consider, in particular their 
locations on the Earth (in terms of northern latitude 5 and eastern longitude a), the 
fiducial masses and the coefficients A and B which appear in the expression of the energy 
resolution, eq. (|59|). 



the neutrinos experience a constant matter potential along trajectories in the mantle and 
a step-like profile mantle-core-mantle along core crossing trajectories. In what follows we 
summarize the analytical results obtained in the two-layers approximation, which cor- 
rectly describe the general features of the Earth matter effects. The comparison with 
exact numerical calculations will be given at the end of this appendix. 

For v e propagating in the mantle only the regeneration factor, f reg , eq. (|23|) , has the form: 

f reg = D{E, 9, 9 n ,p m ) sin 2 (n^J , (60) 

where d = 2R (S cos 9 n is the length of the trajectory in the Earth, i? e ~ 6400 Km is the 
radius of the Earth and p m the matter density in the mantle. The depth D of oscillations 
equals [p~3f : 

D = sin 29 m sin(2# m - 29) , (61) 

with 6 m and l m being the mixing angle and the oscillation length in matter and 9 the 
mixing angle in vacuum. 

Since the mixing is enhanced in matter, 9 m > 9, the depth D is positive as well as the 
whole regeneration factor (|5"0|). From eq. fl5"T| ) one gets: 

sin 2 29 

~ X (x -cos 29Y + sin 2 29 ' (62) 

where 

_ 2EV 
Am 2 



and V is the matter potential. The expression (62) vanishes in the limits of low (x 1) 
and high (s> 1) energies. It reaches the maximum 

Dmax = COS 2 9 , (63) 

at x = 1, which corresponds to the energy: 

E Am 2 

2\f2Gpn e 
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Here Gf is the Fermi constant and n e the electron number density in the medium. 

Thus, the depth of oscillations, D, has a resonant character with E = Er being the 
resonance condition. The width of the resonance is given by the interval E_ 4- E + in 
which D > D max /2. One finds: 

^ = 2 - cos 29 ± J(l - cos 20) (3 -cos 29) , (65) 

which shows that the peak at E ~ is wide for LMA oscillation parameters and gets 
narrower as 9 decreases. Notice that in the limit 9^0 the maximal depth increases, 
D ma x — > 1, but the oscillations disappear due to the decrease of the oscillation phase and 
the vanishing of the resonance width (see eq. (|65|)). 



For antineutrinos, and trajectory in the mantle only, the regeneration factor f reg , eq 
has the same form as in eq. (0), with the oscillation depth 



where we defined 



D = - sin 29 m sin(2# m - 29) , 
sin 2 20 

(x + cos20) 2 + sin 2 20 ' 1 ' 

_ 2E\V\ 
Am 2 

Similarly to the case of v e the depth D, and therefore the regeneration factor f reg , is 
positive, since 9 m < 9\ moreover it has a similar resonant behaviour with maximum 

Dmax = sin 2 9 , (67) 



and the same resonance energy, eq. (|64]) . The borders £L , E + of the resonance interval 
are given by eq. fl65f ) with the replacement cos 29 — > — cos 29. In the limit 9^0 the 
Earth effect disappears due to the vanishing of the oscillation depth flBTp. For maximal 
mixing, cos 29 = 0, we get 



D = D ' 



x 2 + 1 



If the neutrino trajectory crosses both the mantle and the core the analytical treat- 
ment of the regeneration factors, f reg and f reg , is more complicated [p5[| . The interplay 
of oscillations in the mantle and in the core determine irregular oscillations of the factors 
with energy. The depth of oscillations is larger in the region of the energy spectrum close 
to the resonance energies in the mantle and in the core; for SMA oscillation parameters 
parametric effects appear (see fig. |I3p. In contrast to the propagation in the mantle only 
(i.e. in uniform medium), the regeneration factors have negative sign in some intervals of 
energy. 
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Once a realistic density profile of the Earth is considered, the calculation of the re- 



generation factors requires a numerical treatment. The results are presented in fig. |2l 
together with the analytical curves obtained with the two-layers approximation. The fig- 
ure shows the factors f reg and f reg as functions of the neutrino (antineutrino) energy for 
Am 2 = 5 • 10~ 5 eV 2 , sin 2 29 = 0.75 and 6 n = 0°. We used the realistic profile in ref. 
|35| and chose a step-like (two-layers) profile with densities p m = 4.51 and p c = 11.95, 
corresponding to the average densities of the profile of ref. |35 in the mantle and in the 
core along the diameter of the Earth. 

From the figure [H] it follows that the position of the oscillation maxima and minima 
on the energy axis are well reproduced by the step profile. This good approximation of 
the oscillation phase is ascribed to the choice of p m and p c to be the average densities of 
the two layers along the trajectory of the neutrinos. 

In contrast, the depth of oscillations given by the numerical calculation deviates sig- 
nificantly, up to ~ 50%, from the result of the analytical (two-layers) approximation. 
As a general tendency, the depth of oscillations in the realistic density profile appears 
smaller with respect to the case of the two-layers profile. If the density jumps along the 
trajectory are not very large (e.g. for trajectories in the mantle only) this feature can be 
interpreted, qualitatively, by considering the density as smoothly varying along the path 
of the neutrinos. In this case the neutrino conversion occurs adiabatically and the depth 
of oscillations is determined by the matter density at the surface of the Earth. Since the 
surface density is smaller than the average density along the trajectory and the latter in 
turn is smaller than the resonance density pa in the relevant range of energies, a smaller 
depth of oscillations is expected. 

A better approximation of the numerical results can be obtained by using the aver- 
age density in the determination of the oscillation phase and the surface density in the 
determination of the depth of oscillations. 

For SMA parameters the adiabaticity is broken and the two layers model gives a good 
approximation of the numerical results. 



References 

[1] S. P. Mikheev and A. Y. Smirnov, Sov. Phys. JETP 64 (1986) 4. 

[2] S.P. Mikheev and A.Y. Smirnov, Proceedings of the 6th Moriond Workshop Massive 
Neutrinos in Astrophysics and in Particle Physics, 1986, edited by O. Fackler and 
J. Tran Than Van (Editions Frontieres, Gif-sur-Yvette, France, 1986), p. 355. 

[3] A. Y. Smirnov, talk given at the Twentieth International Cosmic Ray Conference, 
Moscow, 1987. 

[4] J. Arafune and M. Fukugita, Phys. Rev. Lett. 59 (1987) 367. 
[5] J. Arafune et al, Phys. Rev. Lett. 59 (1987) 1864. 



30 



[6] H. Minakata et al, Mod. Phys. Lett. A2 (1987) 827. 

[7] P.O. Lagage et al., Phys. Lett. B193 (1987) 127. 

[8] D. Notzold, Phys. Lett. B196 (1987) 315. 

[9] S.P Rosen, Phys. Rev. D37 (1988) 1682. 

[10] T.K. Kuo and J. Pantaleone, Phys. Rev. D37 (1988) 298. 

[11] H. Minakata and H. Nunokawa, Phys. Rev. D38 (1988) 3605. 

[12] L. Wolfenstein, Phys. Lett. B194 (1987) 197. 

[13] A. Y. Smirnov, D. N. Spergel and J. N. Bahcall, Phys. Rev. D49 (1994) 1389 [hep- 
ph/9305204]. 

[14] A. S. Dighe and A. Y. Smirnov, Phys. Rev. D62 (2000) 033007 [hep-ph/9907423] . 

[15] B. Jegerlehner, F. Neubig and G. Raffelt, Phys. Rev. D54 (1996) 1194 [astro- 
ph/9601111]. 

[16] C. Lunardini and A.Y. Smirnov, Phys. Rev. D63 (2001) 073009 [hep-ph/0009356] . 

[17] H. Minakata and H. Nunokawa, Phys. Lett. B504 (2001) 301 [hep-ph/0010240]. 

[18] M. Kachelriess, R. Tomas and J.W.F. Valle, JHEP 01 (2001) 030 [hep-ph/0012134]. 

[19] H. Minakata, Nucl. Phys. Proc. Suppl. 100 (2001) 237 [hep-ph/0101231]. 

[20] K. Takahashi, M. Watanabe and K. Sato, (2000), hep-ph/0012354. 

[21] G. E. Brown, H. A. Bethe and G. Baym, , Nucl. Phys. A 375 (1982) 481. 

[22] H. T. Janka and W. Hillebrandt, Astron. Astroph. Suppl. 78 (1989) 375, Astron. 
Astrophys. 224 (1989) 49. 

[23] S.P. Mikheev and A.Y. Smirnov in ref. 0; for a detailed discussion of regeneration 
of solar neutrinos see e.g. E. Lisi and D. Montanino, Phys. Rev. D56 (1997) 1792 
[hep-ph/9702343]. 

[24] For a review on supernova neutrino experiments see K. Scholberg, Nucl. Phys. Proc. 
Suppl. 91 (2000) 331 [hep-ex/0008044]. 

[25] See W. Fulgione for the LVD collaboration, Nucl. Phys. Proc. Suppl. 70 (1999) 469, 
and references therein. 

[26] See C.J. Virtue for the SNO collaboration, Nucl. Phys. Proc. Suppl. 100 (2001) 326 
[astro-ph/0103324], and references therein. 



31 



[27] The Super-Kamiokande collaboration, Y. Fukuda et al., Phys. Rev. Lett. 81 (1998) 
1158 [hep-ex/9805021]. 

[28] See J.F. Beacom and P. Vogel, Phys. Rev. D60 (1999) 033007 [astro-ph/9811350], 
and references therein. 

[29] A study was presented by the CHOOZ collaboration, M. Apollonio et al., Phys. Rev. 
D61 (2000) 012001 [hep-ex/9906011]. 

[30] The Super-Kamiokande collaboration, S. Fukuda et al., Phys. Rev. Lett. 85 (2000) 
3999 [hep-ex/0009001]. 

[31] The CHOOZ collaboration, M. Apollonio et al, Phys. Lett. B466 (1999) 415 [hep- 
ex/9907037]; Phys. Lett. B420 (1998) 397 [hep-ex/9711002]. 

[32] F. Boehm et al, Phys. Rev. D62 (2000) 072002 [hep-ex/0003022]. 

[33] J.N. Bahcall, P.I. Krastev and A.Y. Smirnov, JHEP 05 (2001) 015 [hep-ph/0103179]. 

[34] A discussion of adiabaticity in the context of supernova neutrinos is given in 
M. Kachelriess and R. Tomas, hep-ph/0104021; see also references therein. 

[35] A. M. Dziewonski and D. L. Anderson, Phys. Earth. Planet. Inter. 25 (1981) 297. 

[36] M.C. Gonzalez-Garcia, C. Pena-Garay and A.Y. Smirnov, Phys. Rev. D63 (2001) 
113004 [hep-ph/0012313]. 

[37] P. Vogel, Phys. Rev. D29 (1984) 1918. 

[38] S.A. Fayans, Sov. J. Nucl. Phys. 42 (1985) 590. 

[39] P. Vogel and J.F. Beacom, Phys. Rev. D60 (1999) 053003 [hep-ph/9903554]; see also 
J.F. Beacom and S.J. Parke, hep-ph/0106128. 

[40] S. Nakamura et al., Phys. Rev. C63 (2001) 034617 [nucl-th/0009012]. 

[41] For review, see e.g. C. Spiering, Annalen Phys. 10 (2001) 131; G.C. Hill, (2001), 
astro-ph/0106064. 

[42] F. Halzen, J.E. Jacobsen and E. Zas, Phys. Rev. D53 (1996) 7359 [astro-ph/95 12080]. 

[43] The AMANDA collaboration, J. Ahrens et al., (2001), astro-ph/0105460. 

[44] A. Bouchta and C. De Los Heros for AMANDA collaboration, private communica- 
tion. 

[45] C. K. Jung, hep-ex/0005046. 

[46] The recent NUSL proposal is available at http://www.sns. ias.edu/~ jnb/. 



32 



[47] A similar criterion was proposed in K. Takahashi et al., (2001), hep-ph/0 105204. 

[48] V. Barger, D. Marfatia and B. P. Wood, Phys. Lett. B 498 (2001) 53 [hep- 
ph/0011251]. 

[49] R. Barbieri and A. Strumia, JHEP 0012 (2000) 016 [hep-ph/0011307]. 
[50] H. Murayama and A. Pierce, hep-ph/0012075. 

[51] For an experimental review see e.g. A. Piepke for the KamLAND colaboration, Nucl. 
Phys. Proc. Suppl. 91 (2001) 99. 

[52] S.G. Wojcicki, Nucl. Phys. Proc. Suppl. 91 (2001) 216. 

[53] A comprehensive review is given in C. Albright et al., (2000), hep-ex/0008064. 
[54] H. Robertson for SNO collaboration, private communication. 
[55] P. Antonioli and W. Fulgione for LVD collaboration, private communication. 
[56] K.S. Hirata et al., Phys. Rev. D38 (1988) 448. 

[57] J.N. Bahcall, P.I. Krastev and E. Lisi, Phys. Rev. C55 (1997) 494 [nucl-ex/96 10010]. 
[58] P. Antonioli et al., Nucl. Instrum. Meth. A309 (1991) 569. 

[59] See the review E.K. Akhmedov, Pramana 54 (2000) 47 [hep-ph/9907435], and refer- 
ences therein. 



33 




Figure 1: The cosines of the nadir angles 9 n of SuperKamiokande, SNO and LVD detectors 
with respect to the supernova as functions of the arrival time of neutrino burst. The 
three panels refer to three different locations of the star in the galactic plane (given by 
the declination angle S s ). We fixed t = as the time at which the star is aligned with the 
Greenwich meridian. „ , 
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Figure 2: Lines of constant flip probability P# in the Am^ tm -tan 2 3 i plane. The solid lines 
refer to E = 50 MeV and correspond, from right to left, to P H = 0.05, 0.1, 0.2, 0.4, 0.8, 0.95. 
The dashed lines correspond to the same values of Pjj with E = 5 MeV. The exclusion 
region from the CHOOZ experiment is shown. 



35 




Figure 3: Lines of constant jump probability in the low density resonance, Pl, in the 
tan 2 Oq-AttIq plane. The solid lines correspond to Pl = 0.05, E = 50 MeV (upper line) 
and E = 5 MeV (lower line). The dashed lines correspond to Pl = 0.5 and, from the 
upper to the lower, E = 40, 20, 10, 5 MeV. The plot of the lines is restricted to the region 
of parameters for which the L-resonance inside the star occurs at densities larger than 
1 g ■ cm -3 (see sect. |2.1| of the text). The allowed regions for the SMA, LMA and LOW 
solutions of the solar neutrino problem are represented. 
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Figure 4: The relative Earth matter effect in v e channel, R, as function of the antineutrino 
energy for LMA oscillation parameters and various values of the nadir angle 9 n . We have 
taken Am| = 5 ■ 10" 5 eV 2 , sin 2 28 Q = 0.75; T- e = 5 MeV, T x = 8 MeV. The figure refers 
to normal mass hierarchy (or inverted hierarchy with Pjj = 1). 
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Figure 5: The same as fig. for 9 n = 60° and various values of AmL 
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Figure 6: The same as fig. ^ for 9 n = 60° and various values of T s 
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Figure 7: The flux factors r and r for various values of the temperatures T e and T g . 
have taken T x = 8 MeV, Ami = 5 • 10" 5 eV 2 , sin 2 26 Q = 0.75 and P H = 1. 
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Figure 8: The same as fig. |] for inverted mass hierarchy and various values of Ph- We 
have taken 9„ = 60°. 
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Figure 9: The relative Earth matter effect in v e channel, R, as function of the neutrino 
energy for LMA oscillation parameters and various values of the nadir angle 9 n . We have 
taken Am% = 5 ■ 10" 5 eV 2 , sin 2 28 Q = 0.75; T e = 3.5 MeV, T x = 8 MeV; P H = 1 (or 
inverted hierarchy). 
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Figure 11: The same as fig. |9] for 6 n = 60° and various values of T e . 
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Figure 12: The same as fig. ^|for 9 n = 60° and various values of Ph- 
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Figure 13: The relative deviations R (solid line) and R (dashed line) as functions of the 
neutrino (antineutrino) energy for LOW oscillation parameters. We have taken Am = 
KT 7 eV 2 , sin 2 20© = 0.9; T e = 3.5 MeV, T- e = 5 MeV, T x = 8 MeV and 6 n = 25°. We 
have also assumed Ph = 1 and Pl = 0. 
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Figure 15: The relative Earth matter effect in v e channel, R, as function of the neutrino 
energy for SMA oscillation parameters and various values of the density profile factor C . 
We have taken Am| = 6 • 10~ 6 eV 2 , sin 2 26 & = 5 • 1(T 3 ; T e = 3.5 MeV, T x = 8 MeV; 
Ph = 1 (or inverted hierarchy); 9 n = 0°. 
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Figure 16: The energy spectra expected at SK, SNO and LVD with (solid lines) and 
without (dotted lines) Earth matter effect, for the same parameters as in figs. ^ and ^| 
and t = 1 hour of fig. |T] a). A distance D = 10 Kpc from the supernova and binding 
energy Eb = 3 ■ 10 53 ergs have been taken. In this specific configuration LVD is not 
shielded by the Earth, thus observing undistorted spectrum. The histogram c) refers to 
the sum of events from z/ e + p — > e + + n and u e + d — > e + + n + n scatterings, while the 
panel d) shows the events from u e + d — >• e + p + p. In a) and b) only the events from 
v e + V — > e + + n are shown. 



49 



t= 1 hour 




1 20 30 40 50 60 70 80 90 

E/MeV 



70 p 

60 r 

50 ~- 

40 r 

30 r 

20 E- 

10 r- 




b LVD 



1 20 30 40 50 60 70 80 90 

E/MeV 



c) SNO V 




1 20 30 



40 50 60 70 80 90 

E/MeV 




1 20 30 40 50 60 70 80 90 

E/MeV 



Figure 17: The same as fig. |TJ] for different values of some parameters: Am^ = 3-10 5 eV 2 
sin 2 26 Q = 0.9; T e = 3 MeV, T g = 4 MeV. 
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Figure 18: The same as fig. ^ for t — 8 hours of fig. p] a). For this configuration SNO is 
unshielded by the Earth. 
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Figure 19: The same as fig. [16| for t — 17 hours of fig. [1] a). For this configuration SK is 
unshielded by the Earth. 
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Figure 20: The ratio of the numbers of low and high-energy events from u e +p — > e + + n 
reaction at SK, SNO and LVD, as a function of the separation energy E s . The bars 



represent 1 a statistical errors. The panels a) and b) refer to the spectra shown in figs. [16 



and \L7\ respectively. We have taken a minimum energy E t h = 20 MeV for the calculation 
of the numbers of low-energy events. 
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Figure 21: The regeneration factors for neutrinos, f reg , and antineutrinos, f reg: calculated 
with step-like (two layers) and realistic profiles. We have taken Am 2 = 5 ■ 10~ 5 eV 2 , 
sin 2 29 = 0.75, 9 n = 0° and the densities p m = 4.51 and p c = 11.95 for the two-layers 
profile. 



54 



